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\—^ Abstract. A dynamical analysis of an effective homogeneous and irrotational 

ff^ Weyssenhoff fluid in general relativity is performed using the 1 + 3 covariant approach 

that enables the dynamics of the fluid to be determined without assuming any 

^ particular form for the space-time metric. The spin contributions to the field equations 

O^ produce a bounce that averts an initial singularity, provided that the spin density 

5—1 exceeds the rate of shear. At later times, when the spin contribution can be 

1 ^V neglected, a Weyssenhoff fluid reduces to a standard cosmological fluid in general 

relativity. Numerical solutions for the time evolution of the generalised scale factor 
^t^ R{t) in spatially-curved models are presented, some of which exhibit eternal oscillatory 

<y-\ behaviour without any singularities. In spatially-flat models, analytical solutions for 

OQ particular values of the equation-of-state parameter are derived. Although the scale 

^'J factor of a Weyssenhoff fluid generically has a positive temporal curvature near a 

bounce (i.e. R{t) > 0), it requires unreasonable flne tuning of the equation-of-state 
L«. parameter to produce a sufficiently extended period of inflation to fit the current 

QQ observational data. 
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1. Introduction 

The Einstein-Cartan (EC) theory of gravity is an extension of Einstein's theory of 
general relativity (GR) that includes the spin properties of matter and their influence 
on the geometrical structure of space-time ([T]; see also [2], [3])- In GR, the energy- 
momentum of the matter content is assumed to be the source of curvature of a 
Riemannian space-time manifold V4. In the EC theory, the spin of the matter has 
been postulated, in addition, to be the source of torsion of a Riemann-Cartan space- 
time manifold f/4 |1] . Weyssenhoff and Raabe [5] were the first to study the behaviour 
of perfect fluids with spin. Obukhov and Korotky extended their work in order to 
build cosmological models based on the EC theory [6] and showed that by assuming the 
Frenkel conditioiJU the theory may be described by an effective fluid in GR where the 
effective stress-energy momentum tensor contains some additional spin terms. 

The aim of this publication is two-fold. First, we wish to investigate the possibility 
that the spin contributions for a Weyssenhoff fluid may avert an initial singularity, as 
first suggested by Trautman [7j. Second, since any realistic cosmological model has to 
include an inflation phase to fit the current observational data, it is also of particular 
interest to see if the spin contributions are able to generate a dynamical model endowed 
with an early inflationary era, as first suggested by Gasperini |8]. Scalars fields can 
generate inflation, but they have not yet been observed. Therefore, it is of interest to 
examine possible alternatives, such as a Weyssenhoff fluid. In contrast to the approaches 
of Trautman [^ and Gasperini [S] , our use of the 1 + 3 covariant formalism enables us 
to determine the dynamics of a Weyssenhoff fluid without assuming any particular form 
for the space-time metric. 

The study of the dynamics of a Weyssenhoff fluid in a 1 + 3 covariant approach 
was initiated by Palle [9]. His work has been revised and extended in our previous 
publication [lOj. The present paper builds on [TO] to extend the work carried out first 
by Trautman [7j in an isotropic space-time, and Kopczynski ^^ and Stewart p2] in an 
anisotropic space-time. It also generalises the analysis of the inflationary behaviour of 
Weyssenhoff fluid models made by Gasperini |8] to anisotropic space-times. 

In our dynamical analysis, we choose to restrict our study to a spatially 
homogeneous and irrotational Weyssenhoff fluid. This particular choice, which implies 
a vanishing vorticity and peculiar acceleration, has been motivated by underlying 
fundamental physical reasons. For a vanishing vorticity, the fluid flow is hypersurface- 
orthogonal, which means that the instantaneous rest spaces defined at each space-time 
point should mesh together to form a set of 3-surfaces in space-time [13]. These 
hypersurfaces, which are surfaces of simultaneity for all the fluid observers, define a 
global cosmic time coordinate determined by the fluid flow. Moreover, by assuming 
that any peculiar acceleration vanishes, the cosmic time is then uniquely defined. It is 
worth mentioning that the absence of vorticity is an involutive property, which means 

I Note that the Frenkel condition arises naturally when performing a rigorous variation of the action. 
It simply means that the spin pseudovector is spacelike in the fluid rest frame. 
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that if it is true initially then it will remain so at later times as shown by Ellis et al 
[H]. Finally, the assumption that there is no vorticity on all scales implies that the fluid 
has no global rotation. This is in line with recent Bayesian MCMC analysis of WMAP 
data performed by Bridges et al. [15] . Their work confirms that a physical Bianchi Vllh 
model, which has a non-vanishing vorticity, is statistically disfavored by the data. 

It is worth pointing out that Szydlowski and Krawiec [16] have considered an 
isotropic and homogeneous cosmological model in which a Weyssenhoff fluid is proposed 
as a potential candidate to describe dark energy at late times. In a subsequent 
publication [17], the authors showed that it is not disfavoured by SNIa data, but it 
may be in conflict with CMB and BBN observational constraints. By contrast, in this 
paper, we consider the full evolutionary history of an, in general, anisotropic universe 
with a Weyssenhoff fluid as its matter source, concentrating in particular on the 'early 
universe' behaviour when the spin terms are significant. Indeed, at late times, when 
the spin contributions can be neglected, the Weyssenhoff fluid reduces to a standard 
cosmological fluid. We thus allow for the presence of a non-zero cosmological constant, 
in accord with current observational constraints. 

In Section [2} we give a concise description of a Weyssenhoff fluid using a 1+3 
covariant approach outlined in Appendix A. The spatial symmetries and macroscopic 
spin averaging procedure are discussed in Section [3j In Section |4| we establish the 
relevant dynamical relations for a homogeneous and irrotational Weyssenhoff fluid. In 
Section [5| we perform a geodesic singularity analysis for such a fluid. In Section [6} 
we analyse the fluid dynamics. The behaviour of the generalised scale factor R{t) of 
such a fluid in a spatially-curved models is discussed in Section^and explicit analytical 
solutions in spatially-flat models are given in Section |8] . For the reader's convenience, 
certain main results obtained in our earlier work [TO] will be repeated in the case of 
a homogeneous and irrotational Weyssenhoff fluid in Section |2] and Section |4j In this 
paper, we use the (+,—,—,—) signature. To express our results in the opposite signature 
used by Ellis [H], the correspondence between physical variables can be found in |10j . 

2. Weyssenhoff fluid description 

2.1. Weyssenhoff fluid phenomenology 

In the EC theory, the effect of the spin density tensor is locally to induce torsion in the 
structure of space-time. In holonomic coordinates, the torsion tensor T^fj_„ is defined as 
the antisymmetric part of the affine connection F'*' , 

which vanishes in GR since the connection is assumed to be symmetric in its two lower 
indices. Note that the tilde denotes an EC geometrical object to differentiate it from 
an effective GR object. 

The Weyssenhoff fluid is a continuous macroscopic medium which is characterised 
on microscopic scales by the spin of the matter fields. The spin density is described by 
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an antisymmetric tensor, 

which is the source of torsion, 

This fluid satisfles the Frenkel condition, which requires the intrinsic spin of a matter 
field to be spacehke in the rest frame of the fiuid, 

Sf^.u" = , (4) 

where u^ is the 4- velocity of the fiuid element. This condition implies an algebraic 
coupling between spin and torsion according to, 

T fiu = i^u Sfj_u , (5) 

and arises naturally from a rigorous variation of the action as shown by ^ . Thus, the 
torsion contributions to the EC field equations are entirely described in terms of the 
spin density. It is also useful to introduce a spin-density scalar defined as, 

S' = IS.^S'^'' > . (6) 

Obukhov and Korotky showed ^ that for a perfect fiuid the EC field equations 
reduce to effective GR Einstein field equations with additional spin terms, and a spin 
field equation. 

The former are found to be, 

where the effective stress energy momentum tensor of the fluid is given by, 

T^. = (ps + PsWu, - PsQ^^u - 2 (/^ + uP'u^) Vp [u(^S,)x] , (8) 

with effective energy density and pressure of the form, 

C2 -lA ^' 

Ps = p — Kb — K A , 

and the physical energy density and pressure satisfy the equation of state, 

p = wp , (10) 

where k = 8ttG, A is the cosmological constant and w the equation of state parameter. 
The spin fleld equation is given by, 

Va (u^S^,) = 2uPu[^\/ix K^pIh) • (11) 
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2.2. Weyssenhoff fluid description in a 1+3 covariant formalism 

The 1 + 3 covariant formalism outlined in Appendix A can now be used to perform a 
more transparent analysis of the Weyssenhoff fluid dynamics. Using a 1 + 3 covariant 
approach in [10], we found that the symmetric stress-energy momentum tensor ([s]) can 
be recast as, 

T^u = PsU^,u^ - Psh^u - 2a^^^S^)x - 2ui^^D^Sy)x , (12) 

where h^^ is the induced metric on the spatial hypersurface, a^^ is the rate of shear 
tensor and D^ is the spatially projected covariant derivative defined in Appendix A. 



Similarly, the spin field equation (11) reduces to, 

S^y + QS^y = 2uPu[^S\p\i,] , (13) 

where B = D^u\ is the expansion rate. 

3. Spatial symmetries and macroscopic spin averaging 

Although much of our following discussion will concern cosmological models that are 
anisotropic, it is of interest to consider the status of a Weyssenhoff fluid as a matter 
source for homogeneous and isotropic models. 

3.1. Spatial symmetries 

To be a suitable candidate for the matter content of such a cosmological model, 
a Weyssenhoff fluid has to be compatible with the Cosmological Principle. In 
mathematical terms, a four- dimensional space-time manifold satisfying this principle 
is foliated by three dimensional spatial hypersurfaces, which are maximally symmetric 
and thus invariant under the action of translations and rotations. 

Although a Weyssenhoff fluid can be expressed as an effective GR fluid, the 
dynamical nature of such a fluid is rooted in the EC theory. Thus, the dynamics of such 
a fluid is determined by the translational and the rotational fields, which are respectively 
the metric g^j_v and the torsion T'^'^j^. The symmetries require the dynamical fields to be 
invariant under the action of an infinitesimal isometry. Hence, the Lie derivatives of the 
dynamical fields have to vanish according to, 

c^g^^u = , (14) 

C-(T\, = . (15) 

where ^^ are the Killing vectors generating the spatial isometrics. A maximally 
symmetric spatial hypersurface admits 6 Killings vectors [18]. The 3 Killing vectors 
^^ generating the infinitesimal translations are related to homogeneity and the 3 Killing 
pseudo-vectors x^ generating the infinitesimal rotations are related to isotropy. They 
satisfy, 

e = h%e , (16) 

X^ = e^f^^D^p^ , (17) 
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where e^''"' is three-dimensional Levi-Civita tensor. 

For a cosmological fluid based on the EC theory, such as a Weyssenhoff fluid, we 
can consider two different forms of the Cosmological Principle: 

(i) the Strong Cosmological Principle (SCP), where the Lie derivatives of the metric 



(14) and of the torsion (15) have to vanish; and 



(ii) the Weak Cosmological Principle (WCP), where only the Lie derivatives of the 



metric ( 14 ) have to vanish and no restriction is imposed on the torsion. 



The translational Killing equation resulting from the symmetries imposed on the 



metric (14) yields. 







1^ 



which is a well-know result obtained in GR. Hence, the WCP is identical to the GR 
Cosmological Principle, which implies that the space-time geometry is described in terms 
of an FRW metric. 



Using the translational Killing equation (18), the rotational Killing equation 



resulting from the symmetries imposed on the metric (15) is found to be 



{D,T\,)e 



(/i-^^TV + h\T^''u + h\T\P)b^^ . 



(19) 



For any maximally symmetric space [18], we can choose respectively a Killing vector 
^^ to vanish at a given point P, and independently, a Killing pseudo- vector x^ to vanish 
at a given point Q according to, 

i^{P) = , (20) 

X'{Q) = . (21) 

Hence, the homogeneity and isotropy can be considered separately. 



By imposing the homogeneity condition (21) on the rotational Killing equation 



(19), the spatial covariant derivative of the torsion tensor has to vanish according to, 

(22) 



DpTx^, = . 

Hence, torsion can only be a function of cosmic time t 
Tx^u = Txf_iu(t) ■ 



(23) 



By imposing the isotropy condition (20) on the rotational Killing equation (19) 



the torsion tensor has to satisfy the constraint, 

h^^h'^^T^, + h^\h''\T^\ + h^rh'^KTx/ = . (24) 

As shown explicitly in a theorem established by Tsamparilis [19] and mentioned 
subsequently by Boehmer [20], the homogeneity (23) and isotropy (24) constraints 



taken together put severe restrictions on the torsion tensor. The only non-vanishing 
components are found to be. 



TxfjLU = hx'^h^ K"'Tlal3-y] = f{t)^\ 



/iiy 1 



rpp 



flp 



u^v^h, 



Ph ^T° 



7/3 



u,,u 



7* 



3"-^ 



T, 



7 5 



(25) 
(26) 
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where f{t) is a scalar function of cosmic time t, p is a fixed index and *T-y is the spatial 
trace of the torsion tensor defined as 

*T, = h^l'T^.p . (27) 

We now discuss the application of this general framework to a Weyssenhoff fluid. 

3.2. Weyssenhoff fluid with macroscopic spin averaging 

The algebraic coupling between the spin density and torsion tensors ^ shows that the 
spin density S^y of a Weyssenhoff fluid can be related to the torsion as, 

^M- = u'^h/K^K-^T^p^ . (28) 



By substituting the non- vanishing components of the torsion (25) and (26) satisfying 
the SCP into the expression for the spin density of a Weyssenhoff fluid (28), it is 
straightforward to show that the spin density tensor has to vanish, 

5^. = . (29) 

Thus, Tsamparilis claims that a Weyssenhoff fluid is incompatible with the SCP 
[19] . This conclusion would hold if all the dynamical contributions of the spin density 



were second rank tensors of the form S^y. However, this is not the case since the 
dynamics contains spin density squared scalar terms. These scalar terms are invariant 
under spatial isometrics like rotations and translations. Hence, they do satisfy the SCP. 
In order for the Weyssenhoff fluid to be compatible with the SCP, the spin density 
tensorial terms have to vanish leaving the scalar terms unaffected. This can be achieved 
by making the reasonable physical assumption that, locally, macroscopic spin averaging 
leads to a vanishing expectation value for the spin density tensor according to, 

{S,,) = . (30) 

However, this macroscopic spin averaging does not lead to a vanishing expectation value 
for the spin density squared scalar since this term is a variance term, 

{S^) = \ {S,.Sn ^ . (31) 

The macroscopic spatial averaging of the spin density was performed in an isotropic 
case by Gasperini [S]. It can be extended to an anisotropic case provided that on small 
macroscopic scales the spin density pseudo- vectors are assumed to be randomly oriented. 

By considering a Weyssenhoff fluid in the absence of any peculiar acceleration 
and by performing a macroscopic spin averaging, we indirectly require the fluid to 
be homogeneous. This follows from the fact that, in this case, the conservation 
law of momentum leads to a vanishing spatial derivative of the pressure and energy 
density. This will be explicitly shown in Section |4.2[ and can also be derived from the 
corresponding dynamical equation for an inhomogeneous Weyssenhoff fluid presented in 
our previous work [TO] . 

Note that even in the absence of a macroscopic spin averaging, the Weyssenhoff 
fluid is still compatible with the WCP, which we discuss further in Section |44 It is 
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worth mentioning that there is no observational evidence so far which would suggest 
that we should impose the SCP even though the mathematical symmetries make such a 
principle mathematically appealing. A true test of whether this principle is applicable 
would be the demonstration of physically observable differences between this case and 
the WCP. 

4. Dynamics of a homogeneous and irrotational Weyssenhoff fluid 

The dynamics of a Weyssenhoff fluid with no peculiar acceleration is entirely determined 
by the symmetric and spin field equations, (JTl), ([l2| and (13) respectively. The former 



can be used to determine the Ricci identities and the energy conservation law. The 
latter simply expresses spin propagation. 

One important consequence of the spatial averaging of the spin density is that the 



stress-energy momentum tensor ( 12 ) reduces to an elegant expression given by 

T^u = PsU^^Uy - pshf,^ , (32) 

where the only spin contributions affecting the dynamics are the negative spin squared 
variance terms entering the definition of the effective energy density and pressure ([9]), 
as expected. These spin squared intrinsic interaction terms S"^ are a key feature that 
distinguishes a Weyssenhoff fluid from a perfect fluid in GR and lead to interesting 
properties we discuss below. 

We have to be careful when performing the macroscopic spin averaging on the 
dynamical equations. The Ricci identities and conservation laws can be entirely 



determined from the stress-energy momentum tensor (32). As we have shown, it 
is perfectly legitimate to perform a macroscopic spin averaging on the stress-energy 
momentum tensor before obtaining explicitly the dynamical equations. However, this is 



not the case for the spin field equations (13). Performing the macroscopic spin averaging 



at this stage would make these field equations vanish. To be consistent, we first have to 
determine the dynamical equations and express them in terms of the spin density scalar 
before performing the spin averaging. 

4-1. Ricci identities 

The Ricci identities can firstly be applied to the whole space-time and secondly to the 
orthogonal 3-space. They yield respectively, 

2V[^V^]np = i?[^^]/MA , (33) 

2Dy^D,^Vp = *i?[Hp^^^A , (34) 

where the spatial vectors v^ are orthogonal to the worldline, i.e. v^u^ = 0, and the 
3-space Riemann tensor *R^upX is related to the Riemann tensor Rf^^px by 

*RnupX = h^^h i,K^ ph xRafS-yS + QfipQuX — Qf^X^up ■ (35) 

The Riemann tensor can be decomposed according to its symmetries as [2T] . 

R'' ux = C^' ux - V^] - ^%^'x} - I'^^V \ ' (36) 
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where C^^ ^^^ is the trace-free Weyl tensor, which, in turn, can be spht into an 'electric' 
and a 'magnetic' part [21] according to. 



E, 



fJ,V 



C, 



fipucr 



P,,'^ 



WU 



H, 



flU 



2^^0-aL' 



aX 



up 



U' 



(37) 
(38) 



The Ricci tensor i?^i, is then simply obtained by substituting the expression (32) for the 
effective stress energy momentum tensor Tf^^^ into the Einstein field equations (JTl), 



-R/, 



V^h 



pv 



(39) 



pv 
''pu - 2 (P'^ + 3Ps) U^U^ - f {ps 

The Riemann tensor Rpp,u\ can now be recast in terms of the Ricci tensor (39), the 
electric (37) and magnetic (38) parts of the Weyl tensor according to the decomposition 
(36) in the following way, 

\k (p, + 3p,) h^\X^Uy^^ - lKp,h^\^h^\^ 



R 



PH 



u\ 



+ AuK^,E^\. - ^h^%,E^\. + 2r/'"^'^n[,if,]. + 27^,^,u^p m^'' . 



(40) 



A] 



It follows from the relation (35) that the Riemann tensor on the spatial 3-space *R' 
becomes, 



pp 



uX 



*R 



pp 



uX 



Khy,h\^Ps - Ah\,E\^ + 2ey,e\^ 



(41) 



The information contained in the Ricci identities ( 33 ) — ( 34 ) can now be extracted by 
projecting them on different hypersurfaces using the decomposition of the corresponding 



Riemann tensors (40) 
pubhcation [TO] . 



(41) and following the same procedure as in our previous 



The Ricci identities applied to the whole space-time yield respectively the 
Raychaudhuri equation and the rate of shear propagation equation. 



= -|0^-2(T^-f (p, + 3p,) , (42) 

(^{pv) = - 1© 0-^1/ - cTi^^ a^)x - E^u . (43) 

The Ricci identities applied to the spatial 3-space express the spatial curvature. 
Their contractions yield the spatial Ricci tensor *TZfj,y and scalar *7l respectively. 



* p 

J^pu 



a 



(pu) 







a 



pu 



-h *TZ 



e^ - 2kPs - 2a' 



(44) 
(45) 



The above expression for the curvature scalar (45) is a generalisation of the Friedmann 
equation. 

One must take particular care when deducing the time evolution of the rate of shear 
from the rate of shear propagation equation (43). This is due to the fact that the rate of 
shear coupling term o-,^au)x and the tidal force term E^^ can not simply be neglected. 
A better route is to deduce the rate of shear evolution equation from the spatial Ricci 
curvature tensor *R^u as shown explicitly by Ellis [22] and outlined below. 

A homogeneous Weyssenhoff fluid satisfies the spatial curvature identity. 



*p li, *'Tp n 

-^pu 3 '^pu '^ ^ • 



(46) 
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Hence, by substituting this identity (46) into the expression for the spatial Ricci tensor 



(44), the propagation equation for the rate of shear is found to be. 



a 



(fiu) 



' flU 



(47) 



This tensorial expression (|47|) can be recast in terms of a scalar relation involving the 

(48) 



rate of shear scalar a according to, 



4-2. Conservation laws 

The effective energy conservation and momentum conservation laws are obtained by 
projecting the conservation equation of the effective stress energy momentum tensor 



(32) 



V^ (i?^. + lg^,n) 



/.v^t;. 







(49) 



respectively along the worldline u'^ and on the orthogonal spatial hypersurface h^j,,^ 
according to, 

-e (ps + Ps) , (50) 

(51) 



Ps 



D^Ps = 



It is worth mentioning that the momentum conservation law (51) expresses the 



homogeneity of the Weyssenhoff fluid. This is due to the fact that according to this 
law, the energy density, the pressure and the spin density of the fluid have to be a 
function of cosmic time only. Hence, the torsion tensor has also to be a function of 



cosmic time only, which is the homogeneity requirement (23). This is only the case for a 



Weyssenhoff fluid with no peculiar acceleration on which a macroscopic spin averaging 



has been performed, as otherwise the momentum conservation law (51) would contain 
additional terms. 

4-3. Spin propagation relation 

The spin conversation law results from twice projecting the antisymmetric field equations 



(13) onto the hypersurface orthogonal to the worldline. 



(52) 



This tensorial expression (52) can be recast in terms of a scalar relation involving the 



spin-density scalar S* in (|6| according to, 

S = -eS. (53) 

This expression implies that the spin density is inversely proportional to the volume 



of the fluid. Note that although the tensorial expression (52) vanishes due to the 



macroscopic spin averaging (30), the scalar expression (53) still applies because it is 



related to the spin variance (31). 
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The effective energy conservation equation (50) can now be recast in terms of tlie 



true (i.e. not effective) energy density and pressure of tfie fluid by substituting tlie spin 



propagation equation (53) 



P = -0(P + P) • (54) 

4.4- Comparision with previous results 

Let us compare our results witli tlie conclusions reached by Lu and Cheng [23] for 
an isotropic Weyssenhoff fluid without any macroscopic spin averaging as shown in 
Appendix A of their publication. 

In an isotropic space-time, the dynamics of a Weyssenhoff fluid, without a 
macroscopic spin averaging, is greatly simplified as we now briefly explain. The 
projection of the effective Einstein field equations ([T]) along the worldline and on the 
orthogonal spatial hyp ersurf aces, yields the following constraint, 

u^/ia^T;, = . (55) 

It arises from the fact that, in an isotropic case, the time-space components of the Ricci 



tensor vanish. From the expression for the stress-energy momentum tensor (12), it is 



clear that the constraint (55) implies a vanishing spin divergence, 

D^S,x = . (56) 

Moreover, the isotropy constraint implies a vanishing rate of shear (i.e. a = 0). Thus, 
in this case, the effective stress energy momentum tensor without the macroscopic 



spin averaging (12) reduces to the elegant expression (32) obtained by performing the 
macroscopic spin averaging. 

Hence, for a Weyssenhoff fluid and isotropic space-time, our results can be compared 
to those of Lu and Cheng [23] • The results of our analysis do not agree with the 
conclusions outlined in [23]. First, they argue that the isotropic Friedmann equation 
implies that the spin density has to be a function of time only, with which we agree. 
Then, they claim that this stands in contradiction with the fact that the spin density 



has also to be a function of space in order to satisfy the projection constraint (55), which 



we dispute. The projection constraint simply implies a vanishing orthogonal projection 



of the spin divergence on the spatial hypersurface (56), which is perfectly compatible 
with the spin density being a function of time only. Hence, contrary to their claim, a 
Weyssenhoff fluid model seems to be perfectly consistent with an isotropic space-time 
(i.e. obeying the WCP), even without spin averaging. 

5. Geodesic singularity analysis 

For a homogeneous and irrotational Weyssenhoff fluid satisfying the macroscopic spin 
averaging condition, the fluid congruence is geodesic. To study the behaviour of such a 
fluid congruence near a singularity, we use the 1-1-3 covariant formalism, which applies 
on local as well as on global scales for a homogeneous fluid model. 
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In order for singularities in the timelike geodesic congruence to occur, the 
Raychaudhuri equation (42) has to satisfy the condition, 
1 



e 



-0<O 



(57) 



near the singularity, as we now explain. First, we recast the singularity condition (57) 
in terms of the inverse expansion rate 6~^ as. 



After integrating with respect to cosmic time t, we find. 



(5J 



(59) 



where 0~^ = 0"^ (t^,) and t = t=K is some arbitrary cosmic time near the singularity. 
Thus, if Q^^ > (0;^^ < 0), the model describes a fluid evolving on a spatially expanding 
(collapsing) hypersurface at t = t,,,. According to the integrated singularity condition 
(59), 0~^ (t) must vanish within a finite past (future) time interval |t — t*| < 3|0~^| 
with respect to t = t^:. Thus a geodesic singularity, defined by 0~^(t) = 0, occurs at 
t = i . 

The homogeneity requirement allows us to define — up to a constant factor — a 
generalised scale factor R according to, 

R 



Q = 3 



R 



(60) 



In a 1 + 3 covariant approach, R is generally a locally defined variable. If the model 
is homogeneous, however, R can be globally defined and interpreted as a cosmological 
scale factor. 

The singularity condition can now be recast in terms of the scale factor R and 
reduces to, 



R 
R 



<0 . 



(61) 



One must also require the scale factor to obey the consistency condition, which requires 
the expansion rate squared to be positively defined at all times according to, 

R 
R 



>0 



(62) 



To determine explicitly these two conditions (61)-(62), the Friedmann (45) and 
Raychaudhuri (42) equations are recast in terms of the scale factor, using respectively 
the expressions for the Ricci (39) and stress-energy-momentum tensor (32) as, 

2 

1 / *'T? 

(r^\ , (63) 




1 / *'JZ 

3 r,„.v + - 



-3 {^^^-^ 



U^U 



-^ + 2a^) 



(64) 
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Using the Friedmann (63) and the Raychaudhuri (64) equations, the consistency (62) 



and singularity (61|) conditions can respectively be explicitly expressed as, 



K 



(^2 + ^+i*^) 



K^^P 



< 1 



S' 



K 



(ff2 _ lA) 1 + 3y; 



-2(^2 



< 



(65) 



(66) 



K,~y 4 

The scaling of the energy density p, the spin density squared S"^ and the rate of 
shear squared a^ can be deduced respectively from the energy conservation law (54), 



the spin propagation equation (53) and the rate of shear propagation equation (48) by 



recasting the expansion rate G in terms of the scale factor i? (60) according to, 

-6 



s' = s' 



a 



cy^ 



p = p 



-3(l+u.) 



(67) 
(68) 
(69) 



Note that the bar corresponds to an arbitrary event (defined by a cosmic time t 
subject only to the condition R j^ 0. 

Furthermore, the spatial Ricci scalar *TZ is the Gaussian curvature of the spatial 
hypersurface, which scales according to, 

*7^ =* n (S\ , (70) 

and the cosmological constant A has by definition no scale dependence, 



A = A ' "^^ 



(71) 



Let us now assume the existence of singularities in the timelike geodesic congruence 
for a homogeneous and irrotational Weyssenhoff fluid. By comparing the scaling 



relations for the spatial Ricci scalar (70) and the cosmological constant (71) with 



those obtained for the spin density squared (67) and the rate of shear squared (68), 
we see that in the limit where the model tends towards a singularity (i.e. i? — > 0), 
the contribution due to curvature and the cosmological constant is negligible. Hence, 
for a Weyssenhoff fluid with a physically reasonable equation-of-state parameter (i.e. 



w < 1), the consistency (65 ) and singularity (66 ) conditions merge into a single condition 
according to. 



«:-V2 



< 1 



K~^p 



(72) 



Moreover, we can recast this condition in terms of the scale dependence R. In the limit 



where the model tends towards a singularity, the condition ( 72 ) becomes, 

-■ill-w) 

< 1 . 



lim 



S^ 



K 



-2^2 



K ^p 



(73) 
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Provided the equation-of-state parameter w < 1, the singularity condition (73) can only 
hold if the rate of shear squared is larger than the spin squared (i.e. a^ > n'^S'^). Hence, 
in the opposite case, where the macroscopic spin density squared of the Weyssenhoff 
fluid is larger than the fluid anisotropies according to, 

K^S^ > a^ , (74) 

there will be no singularity on any scale. This is a generalisation of the result established 
independently for a Bianchi I metric by Kopczynski [11] and Stewart and Hajieck |12j . 
Our singularity analysis is based on the assumption that the Weyssenhoff fluid flow 
lines are geodesies, which implies that the macroscopic fluid (i.e. with spin averaging) 
has to be homogeneous. A key question is whether this still holds in presence of small 
inhomogeneities. According to Ellis [21], the Hawking-Penrose singularity theorems 
apply not only to homogeneous models but also to approximately homogeneous models 
with local pressure inhomogeneities. By analogy, if there is no singularity for geodesies 
fluid flow lines, singularities may still be averted provided the real fluid flow lines can 
be described as small perturbations around geodesies. 



In following sections, we will assume that the spin-shear condition (74) holds, which 
guarantees the absence of singularities for homogeneous models. 



6. Dynamical evolution: general considerations 

For a homogeneous fluid, the Gaussian curvature *TZ depends only on the scale factor 
according to, 

6fc 



*7^ = - 



i?2 



(75) 



where k = { — 1, 0, 1} is the normalised curvature parameter. 

To analyse the dynamics of a homogeneous and irrotational Weyssenhoff fluid, let 



us first recast explicitly the Friedmann (63) and Raychaudhuri (64) equations in terms 



of the physical quantities using the expression for the Gaussian curvature (75) according 
to, 

2 

K [ . 1 / „ ^h M 

(76) 




K 

3 



P 






1 



-A 



(77) 



K \ I 

We will now discuss in more details the geodesic singularities presented in SectzorajS} 
drawing out more fully the geometrical and physical applications. 



6.1. Geometric interpretation oj the solutions 

As outlined above, at stages of the dynamical evolution for which the scale factor R{t) 
is small, a Weyssenhoff fluid with an equation-of-state parameter to < 1 is dominated by 
the spin density and rate of shear contributions. This follows from the scaling properties 
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of the energy (69) of the spin density (67) and of the rate of shear (68). Provided the 
spin-shear condition (74) is satisfied, there can be no singularity (i? -^ 0), because the 



negative sign of the spin squared terms in the RHS of the Friedniann equation ( 76 ) would 



imply the existence of an imaginary rate of expansion, which is physically unacceptable 
(9 G M) as discussed before in Section [s] For physical consistency, the RHS of the 
Friedmann equation has to be positively defined at all times, 



P-kS^ + 



1 



O" 



3k 



+ A >0 



(7J 



which clearly excludes the presence of a singularity provided w < 1. The physical 
interpretation is that — as one goes backwards in cosmic time t from the present epoch 
— the spin contributions to the field equations dominate and produce a bounce, which 
we may take to occur at t = 0, that avoids an initial singularity (i.e. -R(O) > 0). Since 
this model contains no initial singularity, the temporal evolution of the model, governed 



by the Friedmann (76) and Raychaudhuri (77) equations, extends symmetrically to the 



negative part of the time arrow. In order to satisfy the time symmetry requirement and 
avoid a kink in the time evolution of the scale factor i?(t) at t = 0, the expansion rate at 
the bounce has to vanish, -R(O) = 0, and the temporal curvature of the scale factor -R(O) 
has to be finite. Thus, the scale factor R{t) goes through an extremum at the bounce 
i?(0) = i?o dj The energy density at the bounce, po = p(0), is determined by the limit 



where the consistency requirement (|78|) becomes an equality, 

1 



Po 



KDq 



K 



(Tn - 



3A; 
Rl 



A 



(79) 



where 5*0 = S'(O) and ctq = cr(0) denote respectively the spin energy density and the rate 
of shear evaluated at the bounce. Note that this particular choice for the energy density 
(79) at t = has been made in order for the expansion rate to vanish at the bounce. 



This can be shown explicitly by evaluating the Friedmann equation ( 76 ) at the bounce 



using the expression for the energy density ( 79 ) . 



Quantitative expressions or the i?(t)-curve in various cases are derived in SectionM 
below. Before doing so, however, it is worth noting that qualitatively, the general shape 
of the i?(t)-curve for a Weyssenhoff fluid is closely related to the temporal curvature 
of the scale factor R, which is explicitly given by the Raychaudhuri equation (77), and 



also to the range of values for R{t), which is determined by the consistency condition 



(62) on the Friedmann equation (76) 



In this section, let us discuss one particular class of Weyssenhoff fluid models for 
which the cosmological constant A is small (and positively defined), 

< A < Po , (80) 



and the curvature is also small 
< ^ < Po • 

Rq 



^1) 



§ Note that, throughout the paper, a zero subscript denotes the value of a quantity at the bounce (i.e. 
t = 0) and not at the present epoch. 
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The two constraints (80) and (81) on the class of models imply that the sign of the 
temporal curvature of the scale factor depends only on the value of the equation-of-state 
parameter w, which yields three different cases. 



In the first case, where w < — |, the RHS of the Raychaudhuri equation (77) implies 
that the temporal curvature of the scale factor is positively defined at all times, 

R{t)>0 for tG (-00,00). (82) 

The positive sign of R implies that the scale factor is minimal at the bounce and the 
model is perpetually inflating (for t > 0). 



In the second case, where w > 1, by comparing the consistency requirement (78) 



with the Raychaudhuri equation (77), the temporal curvature of the scale factor is found 
to be negatively defined at all times, 

R{t)<0 for te (-00,00). (83) 

Note that for a model with an equation-of-state parameter w > 1, we reach the same 
conclusion as for a fluid with an equation of state parameter w < 1, which is that the 
model has a time-symmetric evolution and bounces at t = 0. The negative sign of R 
implies that the scale factor is maximal at the bounce and is deflating (for t > 0) until 
it eventually collapses. 

In the third case, where — | < if < 1, the symmetric time evolution of the scale 
factor can be split into five parts. Firstly, for a small cosmic time, i.e. |t| < \tf\ — where 
the value of tf depends on the scale parameter w — the sign of the temporal curvature 
of the scale factor is positive. This corresponds to the spin dominated phase. Secondly, 
for a specific cosmic time, i.e. |t| = |t/|, the temporal curvature of the scale factor 
vanishes as the time evolution of the scale factor reaches an inflection point. Then, for 
a larger cosmic time, i.e. |t/| < |t| < \ta\, the temporal curvature of the scale factor has 
the opposite sign until it reaches the second inflection point |t| = \ta\. This corresponds 
to the matter dominated phase. Finally, for large cosmic time, i.e. \t\ > \ta\, the sign of 
the temporal curvature of the scale factor becomes positive again. This corresponds to 
the cosmological constant dominated phase. The behaviour of R{t) in terms of cosmic 
time t is summarised as follows, 

(84) 
(85) 
(86) 
(87) 
(88) 

In the first and second cases, the results obtained for the symmetric time evolution 
of the scale factor are interesting mathematical solutions, but they are inconsistent with 
current cosmological observations. In order to satisfy the current cosmological data, the 
positively defined time evolution of the model has to inflate, i.e. R{t) > 0, at early time 
(t < tf), and produce a sufficient amount of inflation. At later time (t > tf), the energy 



R(t) > 


for 


te{-tf,tf). 


R{t) = 


for 


te{-tf,tf}. 


R{t) < 


for 


te i-ta,-tf)U{tf,ta) 


R{t) = 
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^ ^ \ ^a; ^aj •} 


R(t) > 


for 


t e {-00, -ta) U (ta, 00 
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density of the fluid dominates the dynamics and acts hke a brake on the expansion 
R{t) < 0. 

During the spin-dominated phase, the contribution due to the cosmological constant 



can be safely neglected (80) and the positive temporal curvature of the scale factor (84) 



leads to an inflation phase. The inflatability condition, R{t) > 0, may be deduced from 



the Raychaudhuri equation (77) according to. 



p{l + 3w)-AKS^ + AK-^a^ <0 . (89) 

This inflation phase ends when this inequality is no longer satisfied, which corresponds 
to the inflection point of the temporal evolution of the scale factor, i.e. t = tf. Hence, 
at the end of inflation the density is given by, 

Pf 



(s. 



-2 2 



) 



(90) 



(1 + 3^)^"^ '^ f^ 
The temporal evolution of this model for a positively defined time is characterised 
by a maximal physical energy density p = po coinciding with the start of an inflation 
phase ending when the energy density reaches the density threshold p = pf. At the 
end of inflation, the model enters a matter dominated phase. During this stage, the 
Weyssenhoff fluid model reduces asymptotically to the cosmological solution obtained 
for a perfect fluid in GR in the limit where the cosmic time is sufficiently large t ^ tf, 
which eventually leads to a cosmological constant dominated phase for t ^ ta > tf. 



6.2. Amount of inflation 

The amount of inflation is measured by the number N of e-folds, which is determined 



using the scaling of the energy density (69), the initial (79) and final (90) energy 
densities, and found to be. 



N = \n 



R 



i 



Ro 3(1 +w) 



In 



l + 3w 



.^Sj 



a 



k2 Q2 _ ^2 



(91) 



Using the scaling relations obtained for the spin density squared (67) and the rate 



of shear squared (68), the initial and final values of these quantities are found to be 



related by the number of e-: 


bids according 


to 


S'o = 


-s) 


(Ro^ 
\Rf. 


) = sy , 




^1 = 


-) 


^Ro' 

^Rfy 


' = ay . 





(92) 



(93) 
By recasting the initial values of the spin density squared and rate of shear squared 



in terms of their final values according to (92) and (93) respectively, the expression for 



the number of e-folds (91) reduces to an elegant expression, 

4 



N 



-In 



and is shown in Figure [Tj It worth mentioning that the amount of inflation is 
independent of the rate of shear or the spin density of the fluid. Let us mention that 



Classical big-bounce cosmology 



18 



Bianchi models based on a Weyssenhoff fluid have been studied previously by Lu and 
Cheng [23]. However, the authors did not try to estimate the amount of inflation in 
their analysis. 
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Figure 1. Number of e-folds iV in terms of the equation-of-state parameter w. N{w) 
has a vertical asymptote at w = — ^ 



The only way to have achieve a substantial number of e-folds is by requiring an 
equation of state of the form 

^ - " - (95) 



w 



-3+' 



where < e ^ 1 



which corresponds to no standard fluid and has therefore no acceptable physical basis. 
This conclusion has already been reached by Gasperini [8j in the isotropic case. We 
have showed that the same result still holds in the anisotropic case. 

It is interesting to note that a cosmic string fluid has an equation-of-state parameter 
w = —^. A hybrid Weyssenhoff fluid made for example of fermionic matter cosmic 
strings |25] and matter fields — where the cosmic strings contribution dominates the 



dynamics at the era of interest — has an equation-of-state parameter of the form (95) 



where the value of the fine tuning parameter e depends crucially on the ratio between 
the cosmic string and the matter fields densities. Although such a fluid is a candidate 
to obtain an inflation phase at an early positively defined time (i.e. just after the 
bounce), it does not reduce to the cosmological standard model at later times when the 
spin contribution can be safely neglected. This is due to the fact that the density of 
the cosmic strings contribution pst scales as pst oc R~^. Hence, if the cosmic strings 
contribution dominates the behaviour of the cosmic fluid for an early positively defined 
time, it will do so at all times. 
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However, this problem may potentially be overcome by assuming that, at early 
times, the cosmic strings decay into the matter fields of the standard model leading to 
a reheating phase. It would be worth further investigating this possibility. 

The fine tuning parameter e has a magnitude that is related to the number of e-folds 
according to, 

e ~ e"^^ . (96) 

To obtain, for example, an infiationary phase with A^ = 50 — 70 e-folds — which is a 
characteristic range of values for current parameter estimations — the equation of state 
has to be very fine tuned such that e = 10^*^^ — 10~^^^. It is worth noting that this is 
a similar order of magnitude to the factor 10^^^*^ relating the ratio of the cosmological 
constant predicted by summing the zero point energy of the Standard Model fields up 
to the Planck cutoff to that inferred from cosmological observations, although this is 
almost certainly just a numerical coincidence. 

7. Quantitative dynamical evolution of spatially-curved models 

Our general approach allows one to investigate models with non-zero spatial curvature 
and a cosmological constant. In general, it is not possible to find analytical solutions for 
the time evolution of the scale factor. However, the behaviour of the solutions can be 
analysed by integrating the dynamical equations numerically. The analysis and plots of 
the time evolution of the scale factor in spatially-curved models are presented below. 

7.1. Solutions in presence of a cosmological constant 

The dynamics of a homogeneous and anisotropic Weyssenhoff fiuid in a spatially- 



curved model in presence of a cosmological constant relies on the Fridemann (76) and 



Raychaudhuri (77) equations. Using the scaling relation obtained for the energy density 



(69), for the spin density (67), and for rate of shear (68), the Friedmann (76) and 



Raychaudhuri (77) equations can be recast respectively as. 




-> (^ + ^-' (I) """°' - h' (^» - -^^») (I) " -^ f • (^«' 



R 

where for t = 0, /?o is the scale factor, po the energy density, Sq the spin density and o"o 
the rate of shear. 

For convenience, we introduce six dimensionless parameters defined as, 

R 

r ^ — , (99) 

^t , (100) 
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2 _ ^'^0 
b — , 

npoRl 
A 

KPO ' 



a 



A 



(101) 
(102) 
(103) 
(104) 



which are the scale factor parameter r, the cosmic time parameter r, the rate of 
shear squared parameter 6'^ and the spin density squared parameter s^, the curvature 
parameter a, the cosmological constant parameter A. Note that r and r depend on t, 
whereas S'^, s^, a and A are constant, defined in terms of quantities at the bounce t = 0. 
The consistency condition at the bounce ( 79 ) can be recast in terms of dimensionless 
parameters as, 



6' 



l-a + A 



(105) 



Using (105), the Friedmann (97) and Raychaudhuri (98) equations can also be recast 
respectively in terms of the dimensionless parameters according to, 

.9 1 



( 



3(l-iu) 



a [r 



^r^^l)-l + X(r''-l)) 



3w 



^3(l-«,) ^ ^ 



1-A 



+ 1 



(106) 

, . .^- - .^2 ■-■■■ '^°^' 

where a prime denotes a derivative with respect to the rescaled cosmic time parameter 
r. It worth emphasizing that the dynamics of a homogeneous Weyssenhoff fluid does 
not depend explicitly on the rate of shear parameter squared 6"^. This is due to the fact 
that the rate of shear (68) scales like the spin density (67), and follows explicitly from 
the fact that the 5^— terms cancel after substituting the consistency condition at the 
bounce (105) into the dynamical equations (106) and (107). However, the consistency 
condition implies that the value of the physical quantities at the bounce still depends 
on the corresponding value of rate of shear. 

The physical interpretation of these equations is well known. The Friedmann 
equation corresponds to the conservation law of energy whereas the Raychaudhuri 
equation represents the equation of motion. 

The Friedmann equation ( 106 ) can be recast as follows, 
1/2 ^^ , . a 

where the effective potential is given by 

1 



108) 



f/eff(r) 



2^4 {r'^'-""^ + a-l + \{r'-l)) . (109) 

The parameters present in the Friedmann and Raychaudhuri equations are 
respectively, 

• w: relativistic pressure (SR: continuous parameter). 
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• a: curvature (GR: continuous parameter), 

• —1: spin (EC: discrete parameter), 

• A: cosmological constant. 



From the expression for the effective potential (109), we see that the spin 
contribution has a positive sign, which means that it behaves hke a potential barrier. 
In other words, the spin-spin interaction leads to repulsive centrifugal forces opposing 
the attractive effect of gravity, thus preventing collapse. Note that this is also the case 
for a positive cosmological constant. 

In the absence of relativistic pressure (i.e. w = 0), curvature (i.e. a = 0), spin 
(i.e. the —1 factor vanishes), and cosmological constant (i.e. A = 0) the Friedmann 
and the Raychaudhuri equations reduce respectively to the energy conservation law for 
a particle in a gravitational field with a vanishing total energy (-Etot = 0), and Newton's 
second law of motion. 

The mathematical solutions for the time evolution of scale factor parameter depend 
on the whole real range of the parameters (i.e. w;, a, A G M). But for physical consistency, 
we have to restrict the value of these parameters. Firstly, the Weyssenhoff fluid cannot 
violate causality (i.e. Cg < c), which sets an upper bound on the equation-of-state 
parameter w, 

w<l. (110) 



Secondly, the spin-shear condition (74) and the consistency condition at the bounce 



(105) restrict the range of the cosmological and curvature parameters according to, 

A>a-1. (nil 



In general, it is not possible to flnd analytic solutions for the Friedmann (106) and 



Raychaudhuri (107) equations. However, it is possible to deduce the behaviour of the 
solutions by studying the asymptotic behaviour of the expansion rate parameter r' and 
its derivative r". 

In the limit where r — > 1, the temporal curvature of the scale factor behaves like, 

limr" = -3(|a+|(w;-l)-A) , (112) 

and the expansion rate parameter r' has to vanish, 

limr' = 0, (113) 



to satisfy the consistency condition at the bounce (105). Hence, we flnd three types of 
solutions which depend on the respective value of the parameters: 

(i) lim^^i r" > 0, which implies that the solution r{T) is found within the range 

1 < r < oo , (114) 

provided the parameters w and a satisfy 

X>^a + l{w-l) . (115) 
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(ii) linij,^i r" = 0, which imphes that the solution r{T) is static 
r = 1 , 
when the parameters w and a satisfy 



A 



^a 



+ l{w-l) 
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(116) 
(117) 



A < |q; + ^ (w 



(iii) hnir_^i r" < 0, which imphes that the solution r{T) is found within the range 

< r < 1 , (118) 

provided the parameters w and a satisfy 

1) . (119) 

Moreover, the limit, linir^o ^' = ~oo < 0, clearly does not exist. Hence, the 
solutions always satisfy r > 0, which means that there cannot be any singularity. 
Thus, for a negative temporal curvature r" < 0, the scale factor r reaches a 
minimum value r* found within the range < r* < 1. 

The behaviour of the solutions for the scale factor parameter r{T) is summarised 
in Table [T] below. Explicit numerical solutions in presence of a cosmological constant 
for particular values of the curvature parameter a = { — ^,0, ^} and equation-of-state 
parameter w = { — 1,-^,0,^} are displayed in Figure^- Figure 
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Table 1. Behaviour of the solutions r{w,a,X) 



A 


r 


A>|a+|(«;-1) 


1 < r < oo 


A= la+l{w-l) 


r = 1 


A< la+l{w-l) 


< r* < r < 1 
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Figure 2. [w — —\^a- 
rij) curves for 



Figure 3. (w = — l,a = 0): rij) 
curves for A = { — j^j , 0, 2} 



1.2. Solutions in the absence of a cosmological constant 



In the absence of a cosmological constant, the consistency condition at the bounce (105) 
reduces to, 



s2 - 52 = 1 - a . 



(120) 



Using the consistency condition (120), the Friedniann (106) and Raychaudhuri (107) 
respectively reduce to. 



/2 



( 



3(l-i«) 



a 



{r' - 1) - 1) 



(121) 

, -,- -, , - (122) 

As in the presence of a cosmological constant, the behaviour of the solutions can 
be deduced from the asymptotic behaviour of the expansion rate parameter r' and its 
derivative r". The corresponding results concerning limiting values of r' and r" are 



(l±5£,3a-„,^2(a-l)) 



obtained by setting A = in (112), (115), (117) and (119) . In this simpler case, let 



us now consider the behaviour of the expansion rate parameter r' in the limit where 
r -^ 00 and r —>■ 0. 



(i) w<-l 



„/2 



lim r = 00 > , 



(123) 



which implies that the solutions for the scale factor parameter r{T) diverge 
independently of the value of a. 

(ii) -l<w <1 



lim r'^ 

r— >cxD 



-a > 



(124) 
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Figure 4. (w = ~l,a=i): 
r{T) curves for A = {-|, 0, 2} 
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Figure 9. (w = 0, a = 0): 
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Figure 11. (w=|,a = — ^^ 
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Figure 13. (w=|,a=i): 
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which imphes that the solutions for the scale factor parameter rir^ diverge only 
for a non-closed spatial geometry (i.e. a < 0). Hence for a weakly closed spatial 
geometry (i.e r" > and < a < | (1 — w)), the scale factor parameter oscillates 
between a minimum value r = 1 and a maximum value r* defined by linir^r* r' = 
according to, 

1 < r < r* (125) 

(iii) w < 1, 

„/2 



limr 



-oo < 



(126) 



which does clearly not exist. Hence, the solutions always satisfies r > 0, which 
means that there cannot be any singularity. For a strongly closed spatial geometry 
(i.e r" < and < | (1 — w) < a < 1), the scale factor parameter oscillates 
between a maximum value r = 1 and a minimum value r* defined by limr_>r* r' = 
according to. 



< r* < r < 1 



(127) 



The behaviour of the solutions for the scale factor parameter r{T) is summarised 
in Ta6/e |2] below. Explicit numerical solutions in presence of curvature (i.e. a 7^ 0) for 



w = { — 1, — 35O, 3} are displayed in Figure 14 - Figure 



17 



8. Dynamical evolution of spatially-flat models with zero cosmological 
constant 



In this section, we restrict our study to models with a vanishing spatial curvature and 
cosmological constant (i.e. *TZ = A = 0) and find explicit solutions for the time evolution 
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Table 2. Behaviour of the solutions r{w,a) for A = 



w 


a 


r 


w<-l 


a < 1 


1 < r < oo 




a < 


1 < r < oo 


-l<w<l 


0<«<|(l-w)<l 


l<r <r* 




0<a = l{l-w) <1 


r = 1 




0<|(l-w)<a<l 


< r* < r < 1 



of the scale factor. The reason for choosing this particular class of models is because they 
admit analytical solutions. The dynamics of a homogeneous and anisotropic Weyssenhoff 
fluid in a spatially-fiat model in the absence of a cosmological constant can be solved 
explicitly by determining the asymptotic behaviour of the time evolution of the scale 
factor for particular values of the equation of state parameter. 

It is worth mentioning that exact bouncing solutions for spatially-fiat cosmological 
models based on metric afiine gravity theories (MAG), which include the solutions for 
the class of models based on a Weyssenhoff fluid, have previously been discussed by 
Stachowiak and Szydlowski pU] . 

In the absence of curvature and a cosmological constant, the consistency condition 



at the bounce (105) reduces to, 
s'-6' = l. 



(128) 



Using the consistency condition at the bounce (128) the Friedmann (106) and 



Raychaudhuri ( 107 ) equations can be rewritten in terms of the dimensionless parameters 
according to, 

2 

= r^3(i+-) _ ^-6 ^ (129) 



r 
r 



l + 3w 



'3(l+w) 



+ 2r 



-6 



(130) 



To obtain the explicit time evolution for the scale factor parameter r, the Friedmann 



(129) or Raychaudhuri (129) equations have to be integrated. In order to obtain an 



analytical result, it is easier to integrate the Friedmann equation (129) according to 



r'^dr 



y/r^ii-w) _ I 



d\r\ 



(131) 



The solution of this integral relation depends critically on the value of the 
equation of state parameter w. We will consider six special cases given respectively 

{- 



by w 



■i i 

3' ' 3' 



1,2}, which all admit analytical solutions to (131). The last 
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Figure 14. (w = —X): 

rir) curves for a = {—5,0, j^} 



Figure 15. (w = — |): 
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u: -^ J, 2' ' 4' 2 ' 10 J 



two solutions (i.e. w = 1,2) are physically unacceptable (110) but mathematically 
interesting solutions. 

We first note, however, that in the limit where the model approaches the bounce 
(r — > 1), the asymptotic solution for the scale factor parameter has the quadratic form, 

(132) 



riT)= I l + -{l-wy 

for any equation-of-state parameter w. 

Moreover, in the limit where the model is sufficiently far away from the bounce 
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(r ^ 1), the asymptotic solution for the scale factor parameter is given by, 

riT)= exp(|r|) , (133) 

for an equation-of-state parameter w = —1, and evolves according to, 



/3 \ 3{l+») 

r{T)=\^-{l+w)rj , (134) 

for an equation-of-state parameter w satisfying —l<w<l. Hence, for a positively 
defined cosmic time parameter (r > 0), the asymptotic solutions for the scale factor 
parameter at late times, (r :» 1), have the same time dependence as the solutions found 
within a GR framework. This is due to the fact that the spin contributions can be 
neglected at late times, which implies that the evolution of an effective Weyssenhoff 
fluid asymptotically reduces to a perfect fluid in GR at late times. 

8.1. w = —1 case 

A fluid with an equation-of-state parameter w = —1 behaves like a cosmological 



constant. By solving the integrated Friedmann equation (131) for such an equation- 
of-state parameter, the symmetric evolution of the scale factor parameter with respect 
to the cosmic time parameter is found to be, 

r = (cosh(3r))^/^ . (135) 

For a Weyssenhoff fluid satisfying such an equation-of-state parameter, the symmetric 
temporal curvature of the scale factor parameter 7^{t) is positively defined at all times. 
Hence, for a positively defined cosmic time parameter (r > 0), the model inflates 



perpetually. It starts with a power law inflation phase (132) and tends towards an 



exponentially inflating solution at late times (133). 



8.2. w = —^ case 

A fluid with w = —^ behaves like a macroscopic fluid made of cosmic strings. This 
result has been established by Vilenkin by performing a spatial averaging over a chaotic 
distribution of linear strings made of matter fields [27] . For such an equation-of-state 
parameter, an implicit relation for the symmetric time evolution of the scale factor 
parameter is found according to, 

|r| = iv'?^ + Me(-LF(„,,o,(i),_L)_V2E(a,.ccos(i).-L)) . (136) 

where -F(0, k) and E{(j), k) are the elliptic integral of the first and second kind 
respectively. As in the previous case, the symmetric temporal curvature of the scale 
factor parameter r{T) is positively defined at all times. For a positively defined cosmic 
time parameter (r > 0), the scale factor parameter r{T) tends asymptotically towards 
a constant rate of expansion (i.e. lim i^{T) = 0) in this limiting case. 
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8.3. w = case 

A fluid with w = behaves hke dust. The non-singular behaviour of dust with spin was 
first investigated by Trautman [7] and extended by Kuchowicz [28]. The integrated 



Friedmann equation (131) for an isotropic Weyssenhoff dust can be solved exactly. 



The symmetric evolution of the scale factor parameter with respect to the cosmic time 
parameter is given by, 



-r 



1/3 



(137) 



which agrees with the result established by Trautman. 



8.4. "w = ^ case 

A fluid with w = ^ behaves like radiation. For such an equation-of-state parameter, an 
implicit relation for the symmetric time evolution of the scale factor parameter is found 
according to, 

1 + arccosh (r) ) . (138) 



T 



Vr 



rvr 



As in the anisotropic case, the isotropic solution of the scale factor parameter for a 



relativistic fluid with spin (138) has a clear physical meaning. It is an interpolation 



between two limiting solutions, which describe an inflationary (132) and a radiation 



dominated (134) era respectively. 



8.5. w = 1 case 

A fluid with w = 1 behaves like stiff matter. For such an equation-of-state parameter. 



the derivative of the integrated Friedmann equation (131) with respect to the cosmic 
parameter yields a vanishing rate of expansion, 

r' = . (139) 

The value of the scale factor parameter at the bounce is given by r(0) = 1. Hence, the 
trivial solution for the evolution of the scale parameter with respect to the cosmic time 
parameter is found according to, 

r = 1 . (140) 



8.6. w = 2 case 

Finally, a fluid with w = 2 behaves like ultra stiff matter. A fluid with an equation-of- 
state parameter w > 1 is physically unreasonable given that for such a fluid the speed 
of sound exceeds the speed of light (c^ > c). However, such a solution is mathematically 
interesting because it leads to the presence of singularities. By solving the integrated 
Friedmann equation ( 131 ) for an equation-of-state parameter w = 2, an implicit relation 



for the symmetric time evolution of the scale factor parameter is found according to, 
\t\ = - [ A/r3 (1 - r3) + arctanA/r-3 - 1 j . (141) 
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For a Weyssenhoff fluid with an equation-of-state parameter w = 2, the time- 
symmetric temporal curvature of the scale factor parameter r{T) is negatively defined 



at all times (83). To ensure the continuity of the expansion rate 6 at the bounce, the 



energy density at the bounce po has to satisfy ( 79 ) even if the cosmological solution leads 
to the presence of singularities. As the absolute value of the cosmic time parameter |r| 
increases, the value of the scale factor parameter decreases before eventually collapsing. 



From the implicit dynamical relation ( 141 ), the cosmic time parameter \tc\ at the collapse 



— defined by a vanishing scale factor parameter rdrd) = — is found to be, 

\rc\ = I . (142) 

The collapse of the scale factor i? -^ is equivalent to the divergence of the expansion 
rate — >^ cxd. For an equation of state parameter w = 2, the collapse of the scale factor 
parameter represents a mathematical singularity for the evolution of the scale factor 
parameter with respect to the cosmic time parameter given that the rate of expansion 
diverges at that point, Q{tc) = oo. 

8.7. Graphic solutions 

The cosmological constant {w = —1), cosmic strings {w = — |), dust {w = 0), radiation 
{w = |), stiff matter {w = 1) and ultra stiff matter {w = 2) solutions for the evolution 
of the scale factor parameter with respect the cosmic time parameter r{T) are shown 



in Figure 18 For a positively defined cosmic time parameter (r > 0), the infiection 
point on the graph of r{T) — for the dust and radiation solutions — corresponds to the 
end of infiation. The coordinates of this point are (1.28, 1.41) for the radiation case and 
(1.15, 1.59) for the dust case. 

9. Conclusions 

We have used the 1+3 covariant approach to perform a dynamical analysis of an effective 
homogeneous and irrotational Weyssenhoff fluid. Contrary to the case of a perfect fluid 
in GR, the effective spin contributions to the fluid dynamics act like centrifugal forces 
preventing the formation of singularities for isotropic and anisotropic models satisfying 



the spin-shear constraint (74). The temporal evolution of the models is symmetric with 
respect to t = 0. 

In a cosmological context, the energy density at the bounce state po has 
to be sufficiently dense in order to seed large scale structures from primordial 
quantum fiuctuations. For cosmological parameters which are consistent with current 



cosmological data (80) (81), the temporal curvature of scale factor of a Weyssenhoff 



fluid is positively defined near the bounce (84). However such a fiuid is not a suitable 



candidate for infiation given that the only way to include an infiation phase of about 



50 — 70 e-folds, is by considering a fiuid with a very fine-tuned equation-of-state (95) 
which does not reduce to the standard cosmological fiuid at later times. 
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Figure 18. Symmetric evolution of the scale factor parameter r with respect to 
the cosmic time parameter r for particular values of the equation-of-state parameter 
w = { — 1, — 1,0, |, 1,2} for spatially-flat models with zero cosmological constant. 



It is worth emphasizing that the time evolution of the scale factor of a homogeneous 
and irrotational Weyssenhoff fluid exhibits eternal oscillations, without any singularities. 
By contrast, the corresponding solutions obtained for a perfect fluid in GR are cycloids, 
which do exhibit singularities. Hence, the absence of singularities for a specific range of 
parameters is a genuinely new feature of cosmological models based on a Weyssenhoff 
fluid. 
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Appendix A. 1+3 covariant formalism 

We choose to restrict our scope to a homogeneous and irrotational Weyssenhoff fluid, 
thus implying a vanishing vorticity, cj^i, = 0, and acceleration, a^ = 0. To study 
the dynamics of an such a fluid, we use the 1 + 3 covariant approach which has been 
described in detail in our previous paper [10], and will be summarised below to clarify 
our notation. The approach relies on a 1 + 3 decomposition of geometric quantities with 
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respect to a timelike velocity field u'^ defining an observer according to, 

Vf^u^ = D^u^ = iQh^^ + a^^ = 0^^ , (A.l) 

where 

• hfj,u = QfjLv — u^Uu is the induced metric on the orthogonal instantaneous rest-spaces 
of observers moving with 4- velocity u'^. 

• D^Uu = hffhi/V pUa is the projected covariant derivative of the worldline on the 
orthogonal instantaneous rest-space. 

• O = D^u^ is the scalar describing the volume rate of expansion of the fluid (with 
H = l,Q the Hubble parameter). 

• a^u = Di^^Uy) is the trace-free rate-of-shear tensor describing the rate of distortion 
of the matter flow. 

• Q^u is the symmetric fluid evolution tensor describing the rate of expansion and 
distortion of the fluid. 

It is useful to introduce another scalar quantity, namely the rate of shear scalar defined 

as, 

a' = l^y^ucy^"" > . (A.2) 

Moreover, we define two projected covariant derivatives which are the time projected 
covariant derivative along the worldline (denoted * ) and the orthogonally projected 
covariant derivative (denoted D^). For any general tensor T^''\ , these are respectively 
defined as 

T'^-,. =„^VaT^-,... , (A.3) 

D^T^- .... = h\h^^ ...h%... VeT"-,... . (A.4) 

Furthermore, the dynamics is determined by projected tensors that are orthogonal to 
u^ on every index. The angle brackets are used to denote respectively the orthogonally 
projected symmetric trace-free part (PSTF) of rank-2 tensors and their time derivative 
along the worldline according to, 

rp{f.v) ^ (/^(M^/i^)^ _ Ih^'^hp^) Tf"" (A.5) 

2^W = (/^(M^/i^)^ - Ih^-'hp^) fp'' . (A.6) 

The orthogonal projection of the covariant time derivative of a general tensor T^';^ is 
denoted by, 

{T>^-,J^^h^^...h\...u^Vy^T''-^,„. (A.7) 
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